Unitary equilibrations: probability distribution of the Loschmidt echo 
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Closed quantum systems evolve unitarily and therefore cannot converge in a strong sense to an equilibrium 
state starting out from a generic pure state. Nevertheless for large system size one observes temporal typicality. 
Namely, for the overwhelming majority of the time instants, the statistics of observables is practically indistin- 
guishable from an effective equilibrium one. In this paper we consider the Loschmidt echo (LE) to study this 
sort of unitary equilibration after a quench. We draw several conclusions on general grounds and on the basis of 
an exactly-solvable example of a quasi-free system. In particular we focus on the whole probability distribution 
of observing a given value of the LE after waiting a long time. Depending on the interplay between the initial 
state and the quench Hamiltonian, we find different regimes reflecting different equilibration dynamics. When 
the perturbation is small and the system is away from criticality the probability distribution is Gaussian. How- 
ever close to criticality the distribution function approaches a double peaked, "batman-hood" shaped, universal 
form. 

PACS numbers: 



I. INTRODUCTION 

A sudden change of the parameters governing the evolu- 
tion of a closed quantum many-body system gives typically 
rise to a complex and fascinating dynamics. This so-called 
Hamiltonian quench is attracting an increasing amount of at- 
tention |Q], E St • The reason for such an interest is, at least, 
twofold; in the first place this out of equilibrium phenomenon 
has been recently observed in cold atom systems [5, 6]. Sec- 
ondly, at a more conceptual level, the equilibration dynamics 
of a quenched quantum system plays a role in the very foun- 
dations of statistical mechanics 0, H H [M [U [LI ■ New 
insights can be gained on the fundamental question about the 
emergence of a thermal behavior in closed quantum systems. 

In this paper we study a prototypical dynamical quantity 
for a quantum quench: the Loschmidt echo (LE). This quan- 
tity is defined by the square modulus of the scalar product, of 
the time-evolved, (out of equilibrium) quantum state with the 
initial (equilibrium) one e.g., Hamiltonian ground-state. 

In spite of the simplicity of its definition the Loschmidt 
echo or closely related quantities, convey a great deal of 
information in a variety of physical problems; for example 
X, has been intensively studied in the context of Fermi edge 
sing ularities in x-ray spectra of metals 1X311 . quantum chaos 
OIGJ], decoherence EtnlGjl, and more recently quantum 
criticality [ 19] and out-of-equilibrium fluctuations Il20ll2lll . 

Typically (cfr. figure [TJ the Loschmidt echo rapidly decays 
from its maximum value X - 1 at t = and, after an ini- 
tial transient starts oscillating erraticaly around the same well 
defined value. For finite size systems after a sufficiently long 
time a pattern of collapses and revivals is observed due to the 
almost-periodic nature of the underlying quantum dynamics. 
On the other hand for infinite volume systems the Hamilto- 
nian spectrum generically becomes continuous and an asymp- 
totic value X.OO (coinciding with the average one) is eventually 
reached. 

The main goal here is to investigate the statistical properties 
of the Loschmidt echo seen as a random variable over the ob- 



servation time interval [0, T] . One of the key properties is that 
a small variance, by standard probability theory arguments, 
guarantees that the overwhelming majority of the time X(t) 
sticks very close to its average value J3.E3.II2l- This is the 
sense in which one can speak about equilibration dynamics 
and corresponding "equilibrium properties" of a finite system 
that is evolving unitarily and therefore cannot have attractive 
fixed points. 

We shall show how the features of the probability distribu- 
tion of the Loschmidt echo depend on a rich interplay between 
the initial state and quench Hamiltonian on the one hand and 
the system's size and observation time on the other. In particu- 
lar we will focus on the potential role that the vicinity of quan- 
tum critical points may have on the features of the Loschmidt 
echo probability distribution function H \it\ [l9tl . This latter 
analysis will be mostly carried over by exploiting exact results 
for quasi-free spin chain i.e., the quantum Ising model lfl9ll . 

The paper is organized as follows: in Section|II]we give the 
general setting. Later we introduce a relaxation time Tr and 
discuss the universality content of the X (f) before this time- 
scale. In Section uTCl we define and study other relevant time 
scales, the time T\ for necessary for observing the correct av- 
erage, and revival times where large portion of X (?) are back 
in phase. In section llTDl we give explicit formulas for the mo- 
ments of the LE assuming the non-resonant hypothesis. In 
section [III] we concentrate on a particular example and prove 
all the general results advocated so far for an exactly solvable 
case. Moreover we discover three universal behaviors for the 
whole LE probability distribution function. We draw some 
parallels with another natural quenched observable: the mag- 
netization. Finally section [V] is devoted to conclusions and 
outlook. 



II. GENERAL BEHAVIOR 

Let us start by recalling a few elementary yet crucial facts. 
If H = Yin ^«n» is the system's Hamiltonian (n„'s=spectral 



projections) the closed-system dynamics is described by the 
time-evolution superoperator H t = e~"' H , fiiX) = [H,X]. 
This superoperator is thought here of as a map of the space 
of trace-class operators X into itself Q\X\\\ := trVxtX < 
oo). Closed quantum systems cannot equilibrate in the strong 
sense, as unitary evolutions H, do not have non-trivial i.e., 
non fixed, limit points in the norm topology for t — > oo ||32ll . 
One may then wonder whether a weaker form of convergence 
can be achieved for t — > 00. 

Let us then consider the expectation value of an observ- 
able A(t) := tr(1/ t (po)A) and write the spectral resolution of 
the superoperator H, as a formal sum U, = 2s <?~" £ |£))((£|, 
here £ (|£)))denote the eigenvalues (eigenvector) of 'H, In 
finite dimensions the kernel of ti is spanned by the n„'s 
and gives rise to a time-independent contribution to A(f) 
i.e., Aoo := tr(n„pon„A); the point is now to understand 
whether the remaining components involving the non-trivial 
time-dependent factors exp(-/£f) (£ + 0) admits a limit for 
t — > 00. In finite dimensions the £ are a (finite) discrete set 
of differences of Hamiltonian eigenvalues e.g., £ = E„ - E m , 
and correspondingly A(f) - A„ is a quasi-periodic function: 
the long time limit of A{t) does not exist. On the other hand 
in the infinite dimensional case the spectrum of ti can be 
continuous and, if the function A(£) := (po,£)(£,A) is suf- 
ficiently well behaved, using the Reimann-Lebesgue lemma, 
lim^oo A(£)exp(-/£f) = 0. Therefore in this case 

lim A(t) = A{tj = A 00 (1) 

t— >co 

where A(f) := lim^oo 7 J A(t)dt denotes the time-average 
over an infinite time interval. While this convergence cannot 
be achieved uniformly for the whole set of system's observ- 
ables (in that it would imply strong convergence) it can be 
proven for specific families of A's e.g., local ones 112211 . 

There is a third form of convergence that one can consider 
here: the convergence in probability. In the following we shall 
consider the above defined A(f) as a random variable over the 
real line of t endowed with the uniform measure dt/T with 
T — > 00. Suppose we have a sequence of \Bi\i (think of L 
as the system size), we say that the B^s converge to zero in 
probability if 

Urn Pr{r e R / |B L (r)| > e) = 0, Ve > 0. (2) 

The meaning of this type of convergence should be clear: 
for large L the probability of observing a value of Bi{t) dif- 
ferent from zero is vanishingly small. In other words the 
fractions of f's for which Bi{t) + is going to zero for 
L — > 00. As a matter of fact this is the type of conver- 
gence, with Bi(t) — Ai{t) - Ai co, that has been considered in 
ll 1 OL 11211 . The stochastic convergence (f2]) implies that the prob- 
ability distributions of the random variables Ai i.e., Pz.(a) := 
5(a - Ai(f)) are converging to the one of Aoo := lirriL^co Al.sx, 
i.e., lim^oo Piia) — 5(a — Aca). In this context we say that the 
initial state po "relaxes" or "equilibrates" to p eq if it happens 
thatA^ = Xx{Ap eq ). (13311 . 

A typical strategy to demonstrate this kind of unitary equi- 
libration is to prove (O by showing that i) Bi(t) = ii) 




Figure 1 : Typical behavior of the Loschmidt echo for the Ising model 
in transverse field. All curves refer to a size of L = 100. In the upper 
left panel the relaxation time T R is indicated. Using arguments as 
in section lTtlDK see also notes J23;, 24]) one is able to show that, in 
the quantum Ising model at criticality, the first revival time is exactly 
r rcv = L (upper left panel). 

var(Z?i) goes to zero sufficiently fast for L — > 00. If this is 
the case one can use a basic probability theory result Pr{f e 
R / \Bdt)-B L {t)\ > e) < var(fi L )/e 2 . Since, under the assump- 
tion ii) the RHS of this latter relation can be made arbitrarily 
small and given i), equation (f2]i holds true. 

Notice that yet another way to formulate the kind of con- 
vergence (f2]i is by means of the concept of typicality [8, 9 ]; 
the probability of observing a "non-typical" value i.e. one that 
deviates significantly from the mean one becomes negligible 
in the large L limit. 

A. The Loschmidt echo 

The time dependent quantity we are going to focus on in 
the rest of this paper is the Loschmidt echo (LE): 

£(t) = \me- m \*)f , (3) 

where the state is possibly, but not necessarily, the ground 
state of the Hamiltonian H at a different coupling. In the se- 
quel, statistical averages are always taken with respect to this 
state. 

In the following we will consider X. (0 as a random variable 
with uniformly distributed t > 0. Ideally we are interested not 
only in the first moment but in the whole probability distribu- 
tion function. The probability of X to have value in £2 is given 
by P(£ e Q) = limr^ T- l p{£- 1 (Q) n [0, T]). For those x 
for which the probability density is well defined, it is given by 

P(x) = *(*-£«)= lim 1 J] 1,. 

0<t„<T I dt Khl '\ 

The A:-th moment of this probability distribution is given 
by p k := J x k P(x)dx = -C k (t). Notice that the Loschmidt 
echo can be written as £,(t) = {p^,,e^ fH '(p^)), where: := 
\ifr)(i//\, <H(X) := [H,X] and {X, Y) := ti(X ] 'Y) denotes 
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the Hilbert-Schmidt scalar product. From this it follows 

£"(t) = <p®",e- ,w "" f (p®")), where •H (n) := £" =1 l®^ ® 

® Performing the time average ll34ll one finds p„ = 

(p® n ,'P {n) (p® n )), where !P W projects onto the kernel of ( H i " ) . 
In particular the time average .£ = A*i is given by 

T = ip^^Kp^)) = (P (1) (P{,), r m (p„)> = tr(p e 2 q ) (4) 

where p e? : = !P (1) (p^). From the general discussion in Sect. 
(HHi we know that P W (X) = n„Xn„. The effective equilib- 
rium state p e g = !P (1, (po) is just the po totally dephased in the 
//-eigenbasis. 



B. Short time regime and criticality 

As already pointed out, typically the LE decays from its 
maximum value 1 at t = and, after an initial transient, starts 
oscillating erraticaly around its mean value. In this section we 
will analyze the universality content of this initial transient 
and its dependence on the interplay between the initial state 
\t//) and the evolving Hamiltonian H. 

We start by noticing that the LE Eq. (01 is the square 
modulus of a characteristic function x (0 = (e~" H ) which is 
the Fourier transform of the energy probability distribution: 
X (to) = (6(H — co)). Both x and the LE can be expressed in 
terms of the cumulants of H: 



X (t) = expj]^ (H") c 

n=i 

- (-t 2 ) n 



(5) 



(6) 



where (■),.. stands for the connected average with respect to 
|(/r>. The sums above starts from n — 1 because the zero 
order cumulant is zero: = 0. Since H is a local 

operator, i.e. a sum of local "variables", we can expect in 
some circumstances, the central limit theorem (CLT) to ap- 
ply. More specifically the version of the CLT we are go- 
ing to consider here, is the following. In the thermody- 
namic limit, the probability distribution of the reseated vari- 
able Y = (H — (//)) / -\J(H 2 ) C tends to a Gaussian ( with vari- 
ance 1 and mean zero) . In other words, all but the second 
connected moments ofY tend to zero when the volume goes to 
infinity. 

When the CLT applies, for sufficiently large system sizes, 
the distribution of H will be of the form 



X(oS) = 



1 



exp 



(co - (H)f 



2o- 2 



One can systematically compute corrections to this formula 
and order them as inverse powers of the system size L. Fourier 
transforming back we obtain the characteristic function and 
theLE 



It may seem that this expression for the LE could have be 
obtained right away by keeping only the first term in the ex- 
pansion of the exponential in Eq. (O: 



■fV 2 



X(t) = e"< H >e-2 



£(t) = e- 



(7) 



This is simply a quadratic approximation, that does not rely 
on the CLT. Its validity requires f 2 cr 2 <k 1 that in turn im- 
plies t <s 1 / V(#2) (as will be explained later this means 
roughly t <c Zr''^ 2 where d is the space dimension). From 
figure Q] we see that typically X. (t) decays from 1 and after 
an initial transient, starts oscillating around an average value 
X which will be computed below. Now, when the CLT ap- 
plies, equation (0 can help us to define this transient or re- 
laxation time Tr given by e~ T i>°~ = X- Roughly after this 
time one starts seeing oscillations in X. (?)■ Since in general 
(see below) one has _£ - e ^ L , and when the CLT applies a 2 
scales like the system volume, the relaxation time scales as 

Tr = Y~ In X/o -2 = O The situation is different when 
one considers small variations of the parameters Sh. In this 
limit the average LE is related to the well studied ground state 
fidelity F = |(i/r (1) |i/r (2) )|. More precisely one has X. — F 4 lfl8ll . 
Close to critical points the behavior of the fidelity is dictated 
by the scaling dimension A of the most relevant operator in H 
with respect to the critical state \i[f) [25]. The precise scaling 
is the following: F ~ 1 - const, x Sh 2 L 1{ - d+ ^ h \ where f is 
the dynamical critical exponent. Similarly one can show (see 
below) that <x 2 ~ L 2(d ~ A \ All in all this amounts to saying 
that the relaxation time for small variation 5h around a crit- 
ical point (roughly 5h <k L^ d+ ^ A) ) increases from 0(1) to 
Tr ~ if . In the thermodynamic limit instead, i.e. taking first 
the limit L — > oo, using standard scaling arguments, one can 
show that for 5h small and h close to the critical point h c , the 
relaxation times diverges as T R ~ \h — h c \^ y , v being the cor- 
relation length exponent. In figure [2] we plot the relaxation 
time for a concrete example that will be studied thoroughly 
in section [TTTJ the Ising model in transverse field. There one 
has A = £ = v = 1, and so the singularities observed in fig- 
ure|2]are of simple algebraic type: Tr ~ \h - h c \ around the 
critical points h c = + 1 . 

Let us now turn to discuss when we expect the CLT to work. 
Consider first the case when \i[f) is the ground state of a gapped 
Hamiltonian and the connected energy correlators go to zero 

\x— y\— »oo 

exponentially fast: (H (x) H (y)) — -> {H (x)){H (y)) (expo- 
nential clustering). In this case the connected averages of H 
scale as the volume: (H") c ~ L d . As a consequence the cu- 
mulants of the rescaled variable satisfy (Y n ) c ~ L-(" d - 2d )/ 2 f or 
n > 2, which immediately implies the CLT in the sense given 
above. 

Therefore we can have violation of the CLT only in the gap- 
less case when the state \t//) is critical or when clustering fails. 
Let us then consider a critical state \if/). Connected averages 
have a regular extensive part and a singular part which scales 
according to the most relevant component of H with scaling 
dimension A. When is the ground state of H at a differ- 
ent coupling, A is the scaling dimension of the perturbation 
5H to the critical Hamiltonian. At leading order, we can write 
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Figure 2: (Color online) Relaxation time for the Loschmidt echo 
in the Ising model in transverse field. The critical points are at 
h c = ±1. Clearly we observe divergences at critical points when Sh 
is small. On the line = /z <2) £.(t) = 1 and so there is no relaxation 
or even dynamics. 



f 

Figure 3: (Color online) Rescaled Loschmidt echo JL L (t/ VZ) for the 
Ising model in transverse field. System sizes are, from top to bottom, 
L = 10, 50, 200, 600. The state which defines the average is critical: 
h <r> = 1 while H is at /i (2) = 2.5. The same data collapse feature is 
observed when choosing different h (,) s, although the variance of this 
Gaussian is sensitive to that. 



C. Equilibration and long time behavior 



(H") c ~ A„L d +B n L n(d - &) , and so the rescaled variable satisfies 



(Y n ) c 



A„L d + B n L" id - h) 



(A 2 L d + B 2 L 2( - d -^) 



A)W 2 ' 



If A < d/2 the cumulants of the rescaled variable Y don't go 
to zero but to universal constants J26ll27ll given by 



After having discussed the short time behavior of the LE 
related to the initial transient, let us now turn to its long time 
behavior. 

We first re-write Eq. © in the eigenbasis of H — 
Yj n E„\n){n\ 



-il(E„-E,„) 



(8) 



(Y n ) c 



B„ 

DB/2 - 



In this case the probability distribution of the energy is a, non- 
Gaussian, universal distribution. This kind of universal be- 
havior has been observed for instance in [28] on an example 
where the scaling dimension is A = 1/8. In the opposite sit- 
uation where A > d/2, all the cumulants of Y go to zero ex- 
cept for the first two, and the distribution function approaches 
a Gaussian in the large size limit. In the intermediate case 
A = d/2 the cumulants of Y do not go zero but to a constant 
which is however not universal due to extensive contributions 
coming from the denominator. We recall that in d dimen- 
sional, zero temperature quantum mechanics, operators are 
classified into relevant, irrelevant, and marginal if their scal- 
ing dimension is respectively smaller, larger, or equal to d + £ 
where ( is the dynamical exponent. Hence we see that, even 
in the critical case, we observe deviation from the Gaussian 
behavior only if the perturbation 6H is sufficiently relevant, 
specifically A < d/2. 

To finish let us remind the reader that the CLT also breaks 
down when clustering fails. 

To summarize, when the CLT applies, the LE tends to a 
Gaussian and plotting the function £l yl -^(H 2 ) 2 ) for differ- 
ent sizes L one should observe data collapse (see figure[3]l. 



where p„ = \(iff\n)\ 2 . 

If the spectrum of H is non degenerate the superoperator 
V m acts as a dephasing in the Hamiltonian eigenbasis i.e, 
r m {X) = ^ n {n\X\n)\n){n\. In other words the time average of 
the exponentials in Eq. (O gives simply <5„ m and equation © 
reduces to _£ = 2nP«- As is well known this quantity is the 
purity of an equilibrium, dephased, state: p eq = p n \n){n\. 

Time scales In the preceding section we already defined 
a relevant time scale, the relaxation time Tr which is 0(1) 
off-criticality while T R = O (z/) in the critical case and for 
sufficiently small variations 5h <k L~( rf+ f~ A \ 

In some situations it is useful to consider a finite observa- 
tion time T. We will write £, to indicate the corresponding 
average. It is natural to ask about the interplay between the ob- 
servation time T and the linear size of the system L. In other 
words in general lim L ^oo limj-^co £ L + limj-^co lim^oo £ L . 
Since taking larger system sizes has the effect of sending the 
revival times to infinity and the LE attunes its maximum value 
1 at t = 0, typically the function lim^oo (?) has only one 
large peak at t — whereas £,i (t) has peaks at all the re- 
vival times. Correspondingly we expect lim^co limj-^oo Xl > 
limj-^oo lirn^co £l- This expectation has been confirmed for 
the case of the one dimensional quantum Ising model, see sec- 
tion |m] 

Another question which is relevant in the measurement 
process is how large must the observation time be to effec- 
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T » A 1 , where A 

min' 

trum i.e. A m ; n 



tively measure X? That is, what is the condition to have 
X r ' = X or more in general what is the smallest time T n 

such that one observes [(X)"] 7 " = (X)"? Let us focus on T\. 
Looking at equation ([H) one realizes that it suffices to have 
>, It.-.,-,-. \ j s the smallest gap in the whole spec- 
(E„ - E m ). We can address this ques- 
tion for the class of quasi-free Fermi systems in d spatial di- 
mensions. In this case the energy has the form E n - Y*k n k^k 
where A: is a ^-dimensional quasi-momentum-like label and 
we can assume the one particle energy to be positive. Then 
the gap is given by A min = min A \Y,kPk^k\ with (3 k = 0, ±1. 
By choosing fi k = (-1)* , one obtains a A m ; n which is ex- 
ponentially small in L in those (frequent) cases where A* is 
an analytic function of k J29ll . However in quasi free sys- 
tems the weights p„ decrease exponentially with the num- 
ber of excitations in n. In practice the highest weight is 
given for energy differences between the one and zero par- 
ticle spectra: A <10) = min^ A^ which is a constant of order 1 
in the gapful case, while typically scales as L l for the critical 
case. The next largest amount of spectral weight is attained 

at a gap which is a difference between one particle energies 

A (U) 

= mint g | At - A ? . It will be favorable to have k and 
q nearby in the region where At is flat or almost flat. So we 
get A (u) = min* |A* - K k+dk \ - mint |V*A* ■ 8k\. This gap 
is at least of order of L~ 2 (or at least O (^~ 3 ) if there exists a 
^-vector such that V^A^ = 0). From this discussion we es- 
timate that, at least in quasi-free systems, to have X r > — X 
one must take T\ = 0(1) in the gapful case, while one has 
T\ - O (L) at criticality. If one needs X r ' - X with a larger 
degree of precision than one must choose considerably larger 
time: T x = 0[l 2 ) (or T x = C»(l 3 ) if V k A k = has a solu- 
tion within the allowed set of ^-vectors). Related time scales 
are revival times. We define a revival time to be that particu- 
lar time for which a large portion of spectral weight p n p„, has 
revived. More precisely TVev^peak = 27r, where w pea k is a par- 
ticular frequency E„ - E m such that the weight p„p m is large. 
From the discussion above we expect r rev = (9(1) when H has 
a gap above the ground state, while r rev = O (L) when H is 
critical. These expectations have been confirmed (see figure 
fl} on the hand of a solvable model that will be discussed in 
the next sections. 



D. Moments of the Loschmidt echo 

Having computed the time averaged LE we can now turn 
to higher moments. In doing this one has to distinguish cases 
where n = m from those where E„ = E m in Eq. ([8). So we 
write 



Lit) = £ + X(t), 

X(t) = J]p n p m e- it(E "- E "- } , 



and the n-th moment is given by 



zm n =% n k \r-\x { t)t 

k=0 \ ' 



The computation of the average [X(t)] k can be done assum- 
ing a strong non-resonance condition. With this we mean the 
following. We say H satisfies a ^-non-resonance condition if 
the only way to fulfill Yn=\ Ei, - Ej t — is to match the E/s 
to the Ej's. Strong non resonance is fc-non-resonance for any 
k. Note that this condition cannot be fulfilled when k becomes 
of the order of the Hilbert's space dimension. Now to com- 
pute [X (t)] k draw 2k points in two rows of length k. Imagine 
Ei (Ej) are the points at the left (right). Now draw all pos- 
sible contraction between f s and j's (no contraction among 
Fs or j's since they have the same sing), which are k\, but 
keep only those sets of contractions where there is no hori- 
zontal line. This requirement corresponds to the constraint 
k + ju I = l,...,k. For example for [X(f)] 2 we have only 
one contribution 



so 



For [X(t)] we have two diagrams: 



+ 





Both these diagrams give the same contribution (simply 
swap ;'s with j's) and the result is 



[X(t)f=2 J] p\p\p 



The number of terms in [X(f)] , N (k) is the number of all 
permutations without fixed points and is given by 

TO = £(-r/-'(Mc/!-i). 

7=2 ' ' 

However, among these Af (k) terms, many of them give differ- 
ent contributions. Look for instance at [X (t)] 4 : 

\ 2 



[X(t)f = 3 



Tj p h p l +6 Z pIpIApI 



Correctly one has 3 + 6 = N (4) = 9. 
We collect here the first three moments 

Mi = -C 

P-2 = ^ + YjP\pI 

it +h 



Pi 



= X 3 + 3Xj>^+2 2 p\p\p 
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while the cumulants are 



K\ = Z 

*2 = Yj p l p l 



11*12 



K3 



2 £ PlPlPl 



l\+l2 

h+h-fa+h 



We can notice that each term in [X(t)] has the same form 

— k 

of X. except for a number of non-resonance constraints of the 

1 — k 

form// + i m . Correspondingly [X (t)] < N (k) £ . Using now 

(Yl \ — w 
^ I N(k) - n\ we obtain the simple bound £" < n\£ 

for n>2. This means that 

e XZ <xU)=~^ Z <—^, 
1 ~£A 

and so we obtained a bound on the characteristic function x- 
We see that, when £ — > the probability distribution of £ 
becomes a delta function at zero (the characteristic function 
becomes identically one). The distribution function of the up- 
per bound is 



£ 

Later we will encounter situations where this function gives a 
good approximation to the Loschmidt echo probability distri- 
bution. 



III. ISING MODEL IN TRANSVERSE FIELD 

From now on we will give a detailed description of the 
Loschmidt echo for the case of an exactly solvable model. The 
model we consider is the Ising model in transverse field with 
Hamiltonian 



This model can be mapped to quasi-free fermions and so di- 
agonalized exactly. At zero temperature we distinguish two 
phases: i) An ordered one in the longitudinal direction in 

which (crjOj) — » m , for \h\ < 1, and ii) A paramagnetic 
phase for \h\ > 1. The points \h\ — 1 are critical points where 
the system is described by a conformal invariant field theory 
with central charge c = 1/2. 

The LE is given in this case by Q30I1 (superscript, inserted 
here for clarity, refer to different values of the coupling con- 
stant h) 



-C(t) 



= ]-[(l-sin 2 (^-^)sin 2 (AfV2)) (9) 



where tan (§f) = - sin (k) / + cos (k)) and 
the single particle fermionic dispersion is A^' = 

2 + cos (k)) 2 + sin (A;) 2 . The band minimum 

(maximum) is at E„, = 2 min J|l - /i (2) | , |l + /i (2) |J 

(E M = 2 max {| 1 - h {2) \ , 1 1 + /i (2) |}). Finally, for periodic 
boundary conditions that will be used throughout, the quasi- 
momenta satisfy k„ = n(2n+l)/L,n = 0, 1, . . . ,L/2 — 1. 
Exploiting the fact that H decomposes into a direct sum 
of L/2 blocks 4x4, we are able to compute the complete 
dephased equilibrium state p eq = Pn\n)(n\. The result is 

p eq =Y J P(<*)\a)(a\, (10) 

a-eZ| 

where the multi-index a is a = (a\, . . . , ai), a, = 0, 1, the 
state is \a) = ® k>0 \a k } with 10*} = cos (# (2) /2) |0, 0> A . _ t - 

isin(0< 2 V2)|l,l>*,-* and |0*> = / sin (#< 2 >/2) |0, 0)^ - 
cos [§ (2> /2)\l,l)k,-k- Finally the weights are given by 

r-r tan 2 "' (5^/2) 



k>a 



+ tan 2 (5& k /2) 



(11) 



Using Eq. (TToT > together with (fTTb one can show 
(cfr. ref. rf3llo that the dephased state has the following totally 
factorized form: 



Peq = {£) (a k \0 k }(0 k \ 

k>0 



!»t|l*><l*D 



where a k — {\ + tan 2 (d&k/^)) > an d b k -\ 



a k . 



k>0 



A. Short time regime 

Let us first discuss the short-time, transient regime. Look- 
ing at the function In X. (?) one can readily see that all its n- 
derivatives at t — are the Riemann sums of a summable 
function irrespective of the h^'h being critical. This means 
that all the derivatives of In X. (f) grow linearly with L, and to- 
gether with equation (|6]l implies that the cumulants are linear 
even at criticality, i.e. (h 211 ^ cc L. The same result could have 
been derived by noting that for \h\ + 1 the system is gapful 
and clustering. When is critical, i.e. |/z (1) | = 1, the scal- 
ing dimension of 6H = H (T> - H m = —6h cr\ is one and so 
according to the reasoning in section HlBl the cumulants of H 
grows linearly with L. 

Accordingly the CLT applies. Since all the cumulants, 
including the variance, grow as L, plotting the function 
JLl yl Vl) for different sizes L, one observes data collapse. 
This behavior is illustrated in figure [3] Clearly the plot repro- 
duces a Gaussian with variance 



B. Long time, large sizes and the order of limits 

Consider now a physical situation where an experimenter 
computes X.l- We inserted the labels T and L to stress that 



7 




Since the function £.{t) - e LsW has a limit at infinity, 
its time average is precisely this limit X = e~ Ls(oo) . More 
precisely the distribution function becomes a delta function 
P(x) = s(x-e- Ls{ ° o) ). 

Consider now performing first the time average of Eq. (0. 
According to the discussion in section III CI this requires at 
least observation times as large as T » L (if we are at criti- 
cality). The result in this case is (for the explicit computation 
see the following section) 

2 = exp J] log (l - sin 2 « 2) - 0«) 17) (13) 

fc>0 



ti^ the two 



If now L is large, we can approximate the sum with the in 
tegral: 2 = e ~ L s{ hm < hm ) . Calling 6& k 
functions, g and s (oo) are given by 

g = ~f ln(l-sm 2 (6& k )/2)dk, (14) 



S (oo) 



ln[(l + \cos(6§ k )\)/2]dk. (15) 



We observed that the two averages e~ Lg - obtained by first 
doing the time average and then taking large L - or e~ Ls (°°' - 
obtained by first considering L large and then doing the time 
average - are qualitatively very similar for most values of the 
parameters he". The only region where there is an appreciable 
difference is when hr' and /i <2) correspond to different phases 
(either M 1 -* < 1 and |/! (2, | > 1 or vice- versa). 



Figure 4: (Color online) From top to bottom, mean and variance of 
Loschmidt echo at size L = 100. The region where the variance is 
large shrinks when increasing the system size L. The height of the 
peak however remains constant. 



both size and expectation time are finite, as is required in a 
true experiment. Here we want to study the interplay between 
T and L. Consider first the case where we send L to infinity, 
or more physically L is the largest scale of our system. In this 
situation the spectrum of H is practically continuous, and we 
can write the LE as 



£{t) = exp 



| JT In ( 1 - sin 2 (#f-< 2) ) sin 2 (A»t/2)) dk 



= e -Ls(t) 

In this approximation we sent all the revival times to infin- 
ity and so the function X. (?) is no longer almost periodic, but 
rather tends to a precise limit as t — > oo. We can calculate 
the limit s (oo) and also the first corrections, as t — > oo. The 
procedure is outlined in the Appendix. The result is 

s(t) = s (oo) - j^L cos i^ Em + l n j + (m «_> M ) > (12) 

where .5(00) is the limiting value and A m /M are constants 
which depend on h (l) and are given in the appendix. The re- 
sult d26l i has already been found in [4], here we provide the 
explicit form of the asymptotic value, s (00). 



C. Moments of the Loschmidt echo in presence of degeneracy 

Having the explicit form of the LE Eq. (0 we can com- 
pute its time average and also other moments. Since the Ising 
model is mapped to a free Fermi system on a finite lattice, and 
given the form of the quasi-particle dispersion A*, its spec- 
trum is non-degenerate. In other words the Ising Hamilto- 
nian is 1 -non-degenerate. However, for the same reason, it is 
not £-non-degenerate for k > 2. This means that to compute 
moments higher than the first, we really need to use the ex- 
plicit form Eq. (0 and cannot rely on the results of section 
III CI For the first moment this problem does not arise, and we 
can either use X. = £« p 2 or do the time average of Eq. (0. 
Correctly the results coincide, and they rely on the fact that 
Yjk ( n k ~ m k) At - 0, implies n k = m^, i.e. that the spectrum is 
non-degenerate. The result is 

2=Y\(l-^(^-^)/7). (16) 



k>0 



For later convenience we define a k = sin 2 (flf^ - #[ 2) 
compute higher moments we first rewrite Eq. (0 as 



To 



-C(t) 
X k (t) 



f[(l+X k (t)), 

sin 2 (A k t/2) = J] c*e'^ Ai ' 



S>0 
-a k 



Ok 

2 ' 



3=0, ±1 

Of* 



C ±'= 4 
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We write the n-th power of the LE as 
[£(*)]" = Y](l+Y^(t)), 



A>0 



m=l ' ' y=0,±l,— ,±m 



Now, when computing M products of K*' terms, only the y = 
term will survive after taking the time average. In other 
words: 



/(«) 



..y<») = 



,(n) 
•0,A, 



,(«) 
'0,A- M 



and so we have 



urn =n( i+ *)- 



AM) 



An explicit formula for is 



m=l x ' ft...ftn 



ft. 

ZA=0 



Noting that 

ft,* 
ZA=0 



2 



2«i 



we obtain 



•a - z(t)" 

m=l 



m= i 

For example we have 

s (l) = 

6Q,k 



d?[2t-t 2 - iYl 

' V ' l/=0 



2m 



«A 

' 2 



3 9 

(7i. 4- 

2 



(3) 
•?0,A 



(4) rx ^2^3 4 

£<U = - 2 ^ + 4 a i-4 Q 'A + Y^ Q 'A- 

The variance of the LE is then given by 

A^ = n(l-a,+ ^)-n(l-«*+I^l d7) 

A>0 ' ' A>0 ' ' 

The first moment and the variance are plotted in figure |4] 
One should note that close to the critical points h c — +1 there 
appears a small region Sh where the variance is large. 

Equation dTvT > gives explicitly the variance in a case where 
the non-resonant hypothesis is violated. Since \a^\ < 1 gen- 
erally the variance is given by the difference between two ex- 
ponentially small quantities and so, a fortiori, is exponentially 



small in the system size L. However, looking at figure [4] one 
notes a small region of parameter close to the critical points, 
where the variance is large. As we will see this fact has im- 
portant consequences. 

It is now interesting to compare the result in equation fPTT i 
with what would have been obtained assuming a non-resonant 
condition. Clearly the second moment computed assuming 
non resonance has less terms than the correct one. Since for 
the LE all the contributions are positive, the non-resonant re- 
sult ought to be smaller. In other words, for the variance we 

must have A£ 2 > AX„ r - 

Comparison with the non-resonant result To compute the 
variance assuming non-resonance we use 

& = t + 2j]p 2 p 2 

<<j 

- fi?-Zrf 

i 

Using this formula together with Eq. (fTTti we obtain 

= ^ 2_ n( i_at+ ^) 
n(i-«*+^)-n(i-o»+^)- 



A>0 



We have verified that the inequality A£} > AX.m- holds for 
all values of the coupling constants h (V) and h {2) . However the 

qualitative behavior of AXm- is very similar to the true variance 
AX 1 . 



D. The Loschmidt echo distribution function 

We now turn to consider the whole probability distribution 
of the LE. As we have noted earlier, for any L finite being the 
spectrum discrete, X. ( f) is an almost periodic function. Actu- 
ally X (t) belongs to a smaller class, since it is a trigonometric 
polynomial. In any case, most results we will present are valid 
for the larger class of almost periodic functions. 

We now give the results for the whole LE probability dis- 
tribution function. We have observed three kinds of universal 
behavior emerging in different, well defined regimes, i) Ex- 
ponential behavior where the probability distribution is well 
approximated by # (x) e~ xlL l£, (figure|5]l, ii) Gaussian behav- 
ior, (figure [§}, and iii) A universal double peaked, "batman- 
hood" shaped function, (figure [7]). More precisely we have 
the following scenario: 

• 6h large. In this case, for L moderately large, the dis- 
tribution is approximately exponential. The feature is 
more pronounced when h m , h^ 2 ' are in different phases, 
the limiting case being * ±h ( - x> . 

• Sh small. In this case we have to distinguish two situa- 
tions: 
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- h {y) close to the critical point: 

* L <k - 1 1 ' oc universal batman-hood 
distribution. Note that L <k £ is the so called 
quasi-critical regime. 

* L » | <5/z | 1 , exponential distribution 

- Off critical: 

* L » |<5/ip 2 , exponential distribution 

* Otherwise Gaussian. 

In other words, say that we fixed h> 1 ' in order to have either 
a Gaussian or a double-peaked distribution. We can always 
find an L large enough such that the distribution becomes ex- 
ponential in both cases. However, for the Gaussian case we 
must reach considerably larger sizes L » \8h\~ 2 (compare fig- 
ures [8] and |5). 

We have observed an exponential distribution in the region 
of parameters where the average LE is much smaller than one: 

X « 1. Due to the bound X 2 < 2X , one has AX < X so 
that when £ « 1 even the variance is small. Since the LE 
is supported in [0, 1] and in particular X (?) must be positive 
we expect in the region £ « 1 a distribution with positive 
support, with a large peak very close to zero, and rapidly de- 
caying tail. We have verified that an exponential distribution 

of the form §(x)e^ x ^/X. gives a pretty good approximation 
in the region X « 1. Note in any case, that the exponential 
form is always an approximation. In particular, for x — > the 
true distribution P (x) tends to zero for any value of the param- 
eters. This happens since we have always £(t) > strictly. 
And so generally < Xmin < X < 1. This feature can be 
accounted for by adding a (small) e term to the exponential: 

P (x) oc e -*/2-e/x 

Let us now investigate the conditions under which the first 
moment is small and so we expect an approximately expo- 
nential behavior. Looking at equation ( fl~3T > we see that £« 1 
holds when Lg(h^ ,/z (2 - ) ) » 1 where g is given by Eq. JT4V 

Clearly the function g is zero (its minimum) when h m = h (2 \ 
and is quadratic in the difference on the diagonal. Quite in- 
terestingly the function g is appreciably different from zero 
only when M !) and M 2) correspond to different phases (i.e. ei- 
ther |/z (1) | < 1 and |M 2) | > 1 or vice-versa) so that in these 
cases we observe exponential behavior even for moderately 
small lattices. When h (l> is close to hr\ but away from crit- 
ical points, the function g is quadratic in the difference 6h so 
that X ~ exp (-const, x L5/i 2 ). Hence to have £ « 1 and 
so to observe approximate exponential behavior we obtain the 
relation L » Sh~ 2 . At criticality and for 6h small instead, we 
can use X ~ F 4 where F is the fidelity which scales as (see 
section HlBl and ref. 12511 ') . In the quantum Ising model we are 
considering we have d = ( = v = A= l and so the average 
behaves as X ~ exp (-const, x L 2 5h 2 ^. This means that for 6h 
small around a critical point h c = +1, the condition to have an 
exponential distribution becomes L » 5h' 1 . 

We now turn to consider the origin of the batman-hood 
shaped distribution function. As we have already noticed, the 
LE is a (finite) sum of cosines with given frequencies and 



amplitude. We can imagine a situation where only few fre- 
quencies contribute to the LE. In the limiting case, only two 
non-zero terms. That means that the LE can be approximated 
by 



X (/) = X + A cos (a> A t) + B cos (<o B f) . 



(18) 



where we can assume A, B positive. 

We have devoted some attention to the probability distri- 
bution generated by such a function. If a> A and a> B are ra- 
tionally dependent, the function is periodic and the distribu- 
tion function has square root singularities at all values of X 
where 3,X (0 = 0. However in our case the frequencies oja/b 
are always rationally independent. In this case the vector 
x(t) = (u>At,CL>Bt) wraps around the torus in a uniform way. 
We can then invoke ergodicity and transform the time average 
into a "phase space" average (in this case the phase space is 
x = (x\,Xz)). Hence the probability distribution function is 
given by 

P(£(t)= X )= — J dxi dx 2 6U(M,x 2 )-x)- 
(2n) Jo Jo 

By using Eq. (fT8l this probability density can be written as 



1 r u 



dz 



/max{-lM) ^-^-^(l-a 

(19) 

The integral above can be expressed in terms of elliptic func- 
tions, but we won't need its explicit expression. Typically 
the function ( TT9b is batman-hood shaped function (see fig. |7|, 
with support in |X - \A + B\ , X + \A + B|J and two peaks at 
X — X ± \A — B\. The divergence at the peaks position is of 
logarithmic type, close to the peaks (i.e. x — X ± \A - B\ + e) 
one has 

P Ct) = __^_ + 0( l). 



Note that we never observe the limiting case where B — 
and X (f) becomes a periodic function. This means that even 
a very small spectral weight on B cannot be discarded. On the 
other hand the distribution function ( fT9l ) seems to be quite sta- 
ble against the presence of other oscillating terms with small 
spectral weight. This stability can be seen in fig. [7] where one 
clearly has at least three frequencies with reasonable spectral 
weight, but the probability density is still well approximated 
by a batman-hood. 

Spectral analysis To understand the behavior of P (X = x) 
we do a spectral analysis of X(f) to see which frequencies 
contribute most. In fact, for almost periodic functions there is 
a similar Fourier decomposition as for periodic functions. The 
Fourier expansion is given in this case by Xdisc (&>) = X (0 e l0J '. 
Taking into account Eq. (O Xdisc (&>) can be written as 



Xdisc (W 



(20) 
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Figure 5: (Color online) Approximate exponential behavior. Parame- 
ters are L = 18, h w = 0.3, /i (2) = 1.4. When Sh is large this behavior 
is observed even for moderate sizes (here L = 18), (upper panel). 
The thick line reproduces &(x)e~ x ' J -/X with £, given by equation 
d 1 6b . In the lower panel we plot the Fourier series _£disc (w) given by 
Eq. (2QJ. 



We would like to know which frequencies have the largest 
weight. This is achieved by expanding the product in Eq. (0 

£(t) = l + J]x k (t)+ J] X kl (t)X h (t) + --- 

k>0 k,,k 2 >0 

= £ + Y J X k (t)+ J] X kl (t)X k2 (t) + --- 

k>0 k,,k 2 >Q 

with X k {t) = J] C y h *' = ^ cos (A*0 ■ 

f}=±\ 

Now, since each X k (t) is smaller than 1/2 in modulus, it is 
reasonable to approximate the LE with the first two terms of 
this expansion and we obtain 



X(()-X + JyCO S (A t /). 

k>0 



(21) 



In this approximation we only wrote the zero-frequency con- 
tribution, which corresponds to the mean, and the contribu- 
tion coming from the one particle spectrum. The next term 
has also contributions coming from the two particle spectrum. 
To be more precise, call £ ,(n) the energy of the n-particle spec- 

(first order contribution), while 



trum then E, 



t u oc At, 



7(D 



eP oc a*, 



- b - A fe , and E%> - Eg* oc A*, + A^, (second 

order contribution with less spectral weight). 

Note that we expect Eq. (|2TT > to be approximately valid 
(with a different form for the amplitudes and the frequencies) 



(2) 



oc A, 



Figure 6: (Color online) Gaussian behavior for Sh small but away 
from criticality. Parameters are L = 20, fr 1 ' = 0.1, /i <2) = 0.11. Note 
that the distribution function is an extremely peaked Gaussian (upper 
panel). The thick line is a Gaussian with mean and variance given by 
Eqns. J16t and J17t . In the bottom panel one can notice that many 
frequencies contribute to the LE. The one particle contribution c (w), 
Eq. d22t . is given by the black dots while the red curve gives the true 
spectral decomposition -£ disl . (co) given by Eq. J20b . 



also for non integrable models in which a one-particle approx- 
imation works well. 

Now, if there is a regime where the amplitudes a k /2 are 
highly peaked around few quasi-momenta, in the limiting case 
only two, then the LE can be approximated as in Eq. ( fT8l) and 
we expect a double peaked distribution function. So we are 
led to study the (one-particle) amplitude function 



( \ - ak 



a> e [E m , E M ] . 



(22) 



j=A,- 



Generally c(co) is a bell-shaped function, starting linearly 
from the band minimum E m , reaching a maximum value and 
then decreasing to zero at the band maximum Em- It is not 
difficult to show that [23], when Sh is small and for roughly 
|l -/z (2) | < 10~\ c(a>) starts developing a peak, the width 
of which being proportional to j 1 — h (Z) \. In the limiting case 
M 2) = 1, c(u>) has its maximum at E m = and then de- 
creases monotonically to Em = 2. So, for 6h small and for 
|l - /z (2, | < 10~', c(a>) is a peaked function. In order to have 
few frequencies fall within the peak, and so to have large 
spectral weight on few frequencies, we must additionally have 
L|l - /z (2) | <K 1. This is easily seen analyzing the dispersion 
A k for h (2) close to the critical point 12411 . All in all, the con- 
ditions to have a batman-hood distribution, are 6h small and 
L\l -/i (2) | « 1. The feature is more pronounced when the 
hj are not precisely critical. In fact even though c (a>) is most 
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Figure 7: (Color online) Typical "batman-hood" behavior. Parame- 
ters are L = 40, h w = 0.99, h (1> = 1.01. Upper panel: probability 
distribution (histogram) together with the result of the approximation 
given in Eq. [18] (thick line). The mean £, is taken from Eq. ([9} while 
the parameters A and B are obtained by computing the two largest 
spectral weights (see equation Olt ). Lower panel: the red curve is 
the Fourier series _£ disc (to) (projected spectral density), black curve 
shows the highest coefficient c (to) given by Eq. J221) . together with 
allowed frequencies (black dots). 



peaked when /i (2) = 1 (and the peak is at to — 0), we have to 
remember that the allowed values of a> are a>„ = where 
k„ = n(2n + 1) /L, and the smallest frequency is a>i = A n /i- 
If we perturb h (1) from 1 the peak of c (to) shifts to the right, 
approaching so it is favorable to have h (2) + 1. Not sur- 
prisingly, the conditions to have a batman-hood probability 
density, coincide with having a large variance (see figure [4] 
bottom panel). 

Generally fixing h® and increasing the size L, one eventu- 
ally violates the quasi-critical condition L <K At this stage 
the double-peak feature tends to disappear and the distribution 
approaches an exponential one. This can be clearly seen in fig- 
ure [8] From this figure one can have the impression that the 
"double-peak feature" is a prerequisite of short sizes, since in 
this case one has few frequencies anyway. As we have tried to 
explain instead, this feature survives for larger sizes, provided 
we shrink 5h sufficiently (figure [7]). 

Instead, when 5h is small but fv- 1 ' are far from the critical 
point, then c (at) is not peaked, and many frequencies have a 
large spectral weight (see figure [6]). In this case the distribu- 
tion becomes Gaussian. The emergence of a Gaussian distri- 
bution can be qualitatively understood in the following way. 
First write the LE according to its spectral decomposition 
X. (0 = Yin A n e" a '" where the amplitudes are precisely given 
by A„ = -Cdisc (w„) and are positive. When the frequencies are 




Figure 8: (Color online) Batman-hood distribution approaching an 
exponential one when increasing system size L at fixed h {,) . Parame- 
ters are h m = 0.9, h <2) = 1.2 and chain length are L = 10, 20, 30,40. 



120 




Figure 9: (Color online) Gaussian distribution approaching an expo- 
nential one when increasing system size L at fixed Parameters 
are h m = 0.2, h (2) = 0.6 and chain length are L = 20, 30, 40, 80, 120. 



rationally independent the variables x„ = tu>„ wrap uniformly 
around a large dimensional torus. Then one can consider each 
A„ x e' x " as an independent random variable. The assump- 
tion 5h is small but away from criticality corresponds to 
say that X. (0 can be considered as a sum of many indepen- 
dent random variables, giving rise to a Gaussian distribution. 
When L » |<5/i| 2 the conditions of independence breaks down 
and we recover an approximate exponential behavior. 
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IV. PROBABILITY DISTRIBUTION FUNCTION FOR THE 
MAGNETIZATION 

In the same spirit we can compute the probability distribu- 
tion of a local operator. The first candidate that comes to mind 
is the transverse magnetization. We computed the following 
time dependent observable m (t) — (^f/ (l> \e uH '~' o~e~ l ' Hi ~'\ifr w y 
Using again equation (TTOb one obtains 13611 

m (t) = - V cos (tif) cos (£#fc)+sin sin (6§ k ) cos (/A®) , 

(23) 

where the quasi-momenta range now in the whole Brillouine 
zone: k — n (2n + 1 ) fL, n = 0,1, . . . , L - 1 . Correctly, when 
h m = /z (2) we recover the zero temperature equilibrium result 
<c^> = L- 1 E 4 cos (0*). 

From equation d23l we see that m(t) can be written — 
exactly — as a constant term plus an an oscillating part with 
frequencies given by the single particle spectrum A*. The 
discussion on characteristic times becomes simplified as all 

(2) 

time scales are uniquely determined by A^ . For example 

the time T\ necessary to observe the correct mean: m T * = m 
must simply satisfy T\ » gap -1 which means T\ » L in the 

quasi-critical regime |M 2) — l| » L, while it suffices to have 
T\ » 0(1) away from criticality. Given equation (1231 it is 
not difficult to compute the mean and the variance, which are 
given by 

m = - V cos (j?^) cos (5^) (24) 
L k 

Am 2 = J_ ^ sin 2 (#f ') sin 2 (5§ k ) . (25) 

Some comments are in order here. First fixing h <l \ /i (2) the 
variance (figure |T0T> goes to zero as Zr 1 and not exponentially 
fast as was the case for the LE. Second, the spectral weight 
associated to the frequency A^ is sin (#< 2) ) sin (6$ k ). We ob- 
served that there are always many frequencies with large spec- 
tral weight. In other words the spectral weight function is 
never peaked. 

These comments suggest us to expect a Gaussian behav- 
ior for the probability distribution function of the magnetiza- 
tion irrespective of the parameters approaching critical values. 
This has indeed been observed (fjg.fTlTi. 

V. CONCLUSIONS 

The unitary character of the dynamics of a closed quantum 
system implies that whatever relevant notion of equilibration 
one might have has to be a subtle one. In this paper we inves- 
tigated the unitary equilibration of a quantum system after a 
sudden change of its Hamiltonian parameter. To this aim we 
used a prototypical time-dependent quantity: the Loschmidt 
echo. We established how the global features of _£ depend 
on the physical properties of the initial state preparation and 




Figure 10: (Color online) Variance of the magnetization as given 
by Eq. dlOt for L = 80. A signature of criticality are the cusps at 
= ±1. 
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Figure 1 1 : (Color online) The probability distribution for the mag- 
netization has only Gaussian behavior. Here parameters are L = 
40, h m = 0.9, h {2) = 1.01. The continuous line is a Gaussian with 
mean and variance given by Eqns. d24l > and l |25| l. 

on those of the quench Hamiltonian. The central object of 
our analysis is given by the long time probability distribution 

for £: P(x) = SU(t) - x) := lim^ T~ l £ 6(£(t) - x)dt. 
Broadly speaking concentration phenomena for P correspond 
to quantum equilibration. 

Here below for the reader's sake we summarize the main 
findings of the paper 

• Resorting to a cumulant expansion we characterized the 
"short" time behavior of -C(t). Different regimes can be 
identified depending on the most relevant scaling di- 
mension of the quench Hamiltonian. When the central 
limit theorem (CLT) applies one has a Gaussian decay 
over a time scale 0(1) for gapped systems. For the criti- 
cal case the time-scale becomes O (z/) in a small region 

\5h\ <s L- (d+( ~ A \ At critical points the CLT can be vio- 
lated and X.(t) takes a universal non-Gaussian form for 
sufficiently relevant perturbations. 

• We discussed the general structure of the higher mo- 
menta of P i.e., ji k :- X. k (t) = J x k P(x)dx using the 
so-called non-resonant hypothesis. We showed that all 
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the Hk are bounded by those corresponding to a Poisso- 
nian distribution i.e., P(x) = -&(x)exp{-x /£)/£. 



Using exact results for the quantum Ising chain we rig- 
orously analyzed the interplay between the chain length 
L and the averaging time T. In particular we showed 
how the limit lim^oo and the thermodynamical one 
i.e., lim^oo do not commute. While in finite sys- 
tems the .£(0 is an almost-periodic function, in the 
thermodynamical limit Xco := lim,^ co X(f) exists and 
P(x) — > 6(x - Xco). We explicitly computed Xoo and the 
way it is asymptotically approached for large t . We gave 
a general closed form for the exact fi^'s and compared 
with that obtained with the non-resonant hypothesis 



For the quantum Ising chain we numerically investi- 
gated P(x). We identified three universal regimes a) An 
exponential one (P is Poissonian) when L is the largest 
scale of the system b) A Gaussian one for intermedi- 
ate L and initial state and quench parameters close and 
off-critical c) A "Batman-hood" shape for P when the 
parameters are close to each other and close to crit- 
icality. This result holds in the quasi-critical region 
L\h & - 1|«C 1, (2 = 1,2). 



Finally, for the sake of the comparison of the Loschmidt 
echo with a prototypical observable, we computed the 
time-dependent magnetization after the quench and 
studied its long-time statistics. In this case only a Gaus- 
sian regime appears to be reachable. 



We have shown that the Loschmidt echo encodes sophisticated 
information about the quantum equilibration dynamics. For fi- 
nite system vastly different time-scales arise: short time relax- 
ation is intertwined with a complex a pattern of collapses and 
revivals and eventually Poincare recurrences. Unveiling how 
these phenomena depend on spectral properties of the under- 
lying Hamiltonians, is one of the key challenges in the way 
to understand emergent thermal behavior in closed quantum 
systems. 
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Figure 12: Typical behavior of the function s(t) (black line) in the 
thermodynamic limit. The red line is the approximation given by 
Eq. {26]». Parameters are h (t) = 1.3, h {2) = 2. 



Appendix 

Asymptotic of s (?) 

Here we want to compute the asymptotic of s (f) for t — > oo, 
that is the integral: 

S W = ~Yn f ^ t 1 " Si " 2 W° _ **°) Si " 2 ( A ^ /2 )) dk 

(2) 

We can go to energy integration setting A k — to 
1 C Em 

s(t) = -— J ln|l -<?(w)sin 2 (wf/2)Jp(w)fiw. 

Where E m = 2 min {|l + h (2) \ , |l - h (2) \) and E M = 
2 max {1 1 +/z (2) |,|l -h (2) \]. To be explicit: 



a{oS) 



2a> 



^co 2 -El)(E 2 M -u 2 ) 

(^-j%)(j%-^)(ftg>-ftt'>) 2 
4 (M 2 >) 2 [4 (M 2 > - MD) (1 - A W>) + M«w 2 ] a?- ' 



Note that a («) is zero at the band's edge, positive otherwise 
(and smaller than 1 in modulus). Instead p (a>) has square root 
(van Hove) singularities at the band edges as a result of the 
quadratic dispersion at those points (when M 2) + 1). When 
/z (2) = 1 the dispersion is linear at the bottom of the band but 
still quadratic at the upper band edge, hence in this case only 
the square root singularity at the upper band edge survives. 

Then expand the logarithm into an infinite series. Using the 
Riemann-Lebesgue lemma we can show that, 



lim f 



f(w) [sm(wt)] 2k dw = 2 



■:/: I 2£ 
k 



)! 



f (w) daj, 



provided that / is summable. The resulting series can be 
summed 



-Z 



-2k 



2k 



= 2 In 



(i+ VT^ 



if Ixl < 1. 
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So finally 

1 r 

lim s (?) = — 

n Je„, 

1 r 

n Jo 



£m [ 1 + VI -a(oj) 
In 



p (a>) da> 



In 



1 + VI -a(k) 



dk. 



To compute the first correction to the limit note that a (w) 
smoothen the singularity at the band edge, so that for the lead- 
ing correction we need only k — 1 in the expansion of the log- 
arithm. The evaluation of the oscillating integral is done with 
a saddle point technique. The result is 



S (t) S (oo) • 



i r 

— I a (co) p (<d) cos (tot) du> 
4jT Je,„ 



s (oo) - 4^ cos UE m + ln] + lm <•-> A It. ( 26 ) 



with constants given by (we assumed here h {2) > so that the 



band minimum is E m = 2 1 - A (2) |) 

1 (a« - M 2 >) 2 



Am = - 



16 V* (1 - AC 1 )) 2 (AC 2 )) 3/2 J|l-AC 2 )| 
1 - A (2) ) 2 



16 V* (1 + AC 1 )) 2 (A( 2 )) 3/2 JT+A^j 

The result d26b should be the same as the square modulus of 
Eq. (12) in |4|. However in J3l there appears only one fre- 
quency, corresponding to the lowest band edge. The discrep- 
ancy probably arises from a continuum approximation which 
discards the effect of the van Hove singularity present at the 
upper band edge. As we have seen both terms give similar 
contributions. In particular, even at criticality, the van Hove 
singularity at the upper band edge survives. 

In figure[T2]one can appreciate the validity of the approxi- 
mation I 
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